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Abstract 

We report on the first detection of maser emission in the J=ll-10, J=14-13 and J=15-14 transitions 
of the v—Q vibrational state of SiS toward the C-rich star IRC+10216. These masers seem to be 
produced in the very inhomogeneous region between the star and the inner dust formation zone, 
placed at ~5-7 R*, with expansion velocities below 10 km s~^. We interpret the pumping mechanism 
as due to overlaps between v—1-0 ro-vibrational lines of SiS and mid-IR lines of C2H2, HCN and 
their ^■^C isotopologues. The large number of overlaps found suggests the existence of strong masers 
for high-J v—Q and v—1 SiS transitions, located in the submillimeter range. In addition, it could be 
possible to find several rotational lines of the SiS isotopologues displaying maser emission. 
Subject headings: circumstellar matter — masers — stars: AGB and post-AGB — stars: carbon 



1. INTRODUCTION 

The detection of strong maser emission at the frequen- 
cies of pure rotational transitions of some molecules is a 
common pheno menon in circumstellar env elopes (CSE's) 
of evolved stars (|Elitzuilll99i lGravlll999f ). The maser is 
usually produced in a small region of the envelope and 
sometimes provides valuable information on the physical 
conditions of the emitting region. 

Due to the different chemistry, masers are produced 
by different molecules in O- and C-rich stars. In 0-rich 
stars, SiO exhibits strong maser emission in different 
rotational transiti ons within severa l vibra t ional states, 
from v=\ to 4 (jCernicharo et al.l 119931 : iPardo et al.l 
[1991 . These masers are formed in a region of the 
CSE very close to the st ellar surface and seem to be 
driven by NIR radiation (jPardo et al.ll2004D . In C-rich 
stars, although SiO i s present with sim ilar abundances 
than in 0-rich stars (jSchoier et al.l[2006h . no SiO maser 
has been detected. The explanation could be that SiO 
is formed at ~3-5 R*, where the angular dilution of 
the star is high and the density and temperature lower 
than in the regions whcrc_SiO maser s are produced in 
0-rich stars figundez fc Cernicharol I2006D . In C-rich 
stars only HCN shows strong maser emission in sev- 
eral pure r otational lines within vibrational states from 
V2=l to 4 (|Lucas fc Cernicharol ri989t iSchilke fc MentenI 
[200I . These masers must be formed in the innermost 
regions of the CSE. SiS has been previously found to 
show weak mas er emission in the J= 1-0 v=Q transition 
in IRC-hl0216 (iHenkel et al.lfT98l . 

In this letter, we report on the first detection of maser 
emission from the J=ll-10, 14-13 and 15-14 transitions 
in the v—0 vibrational state of SiS (hereafter Mi , M2 , and 
M3) observed toward the C-rich star IRC-f 10216. We 
have also obtained observations of v—\ rotational lines 
which exhibit thermal emission. We propose that over- 
laps of w=l-0 ro-vibrational transitions of SiS with mid- 
IR lines of C2H2 and HCN could provide the pumping 
mechanism for these masers as well as higher- J v—Q SiS 
masers in the submillimeter range. This discovery is in- 
teresting because this species could play in C-rich stars a 
role similar to that of SiO in 0-rich stars: the energy level 
pattern of both molecules is similar and it is also formed 



close to the star, as chemica l equilibrium and interfer- 
ometric observations imply (jBieging fc Nguven-O-Rieiil 
[T989h . 

2. OBSERVATIONS 

The observations of the v=0 J=6-5 and J=8-7 to 
J=15-14 transitions of SiS (see left panels, Fig. [T|) were 
carried out on 2004 June 19**^ with the IRAM 30 m ra- 
dio telescope. Four SIS receivers operating at 1, 1.3, 
2, and 3 mm were used simultaneously. System tem- 
peratures were in the range 120-225 K for the 1, 1.3 
and 2 mm receivers and 200-600 K for the 1 mm re- 
ceivers. Atmospheric opacities ranged between 0.08 at 
108 GHz to 0.28 at 267 GHz. For the J=10-9 SiS line 
at 181.5 GHz, the system temperature was significantly 
higher, ~10* K, due to proximity to the atmospheric 
water line at 183.3 GHz. The intensity scale was cal- 
ibrated using two absorbers at different temperatures 
according to the At mosphere Trans mission Model ATM 
(|Cernicharo 1985; Pardo et al.ll200ll ). Pointing and focus 
were regularly checked on the nearby quasar OJ 287. The 
observations were made in wobbling mode, with 180" off- 
set and the secondary nutating at a rate of 0.5 Hz. The 
backends were a filter bank with 256 MHz bandwidth and 
1 MHz resolution and an autocorrelator with 80 kHz res- 
olution (Au~0.08-0.2 km s"^). The lines of SiS in the 
v=l state (see right panel, Fig. [T|) were also observed 
with the IRAM 30 m telescope between 1989 and 2001. 

Additionally we show some spectra corresponding to 
ro-vibrational lines of SiS (see center panels in Fig. [1]) 
from a mid-IR high resolution spectrum (11-14 /im) 
obtained in 2002 December with the TEXES spec- 
trometer mounted on the 3 m tele scope IRTF (ob- 
servational procedures desc ribed in iLacv et all 120021 : 
iFonfrfa Expdsito et al.ll2006f) . 

3. RESULTS AND DISCUSSION 

The v—1 rotational lines show single cusped profiles 
and their relative intensities indicate that the rotational 
levels are thermally populated. The linewidths corre- 
spond to velocities of 9-11 km s~^, lower than the ter - 
minal velocity ~14.5 km s~^ (|Cernicharo et al1l2000f) . 
thus the emission arises from the innermost region of the 
CSE, between the photosphere and the inner dust for- 
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Fig. 1. — Millimeter and mid-IR observations towards IRC+10216. - The pure rotational transitions in the v=0 state (left inserts) and 
in the v=l state (right inserts) are plotted as a function of the velocity with respect to the star (vijSfl=— 26.5 km s~^). The velocity 
structure of the three observed masers Mi, M2 and M3 (plotted as red histograms) is shown in more detail in the middle-bottom insert, 
where the upper lines have been increasingly shifted by 10 K. — The middle inserts present some mid-IR spectra at the frequencies of 
several ro- vibrational v=l-0 SiS lines, whose labels and frequencies are also shown. The zero of the velocity scale is set to the SiS frequency 
corrected for the LSR. The observed spectrum is plotted as a black histogram . R ed lines correspond to a radiative transfer model which 
only considers C2H2, HCN and their isotopologues (Fonfn'a Exposito et al. 20061) . The vertical lines indicate the position of maximum 
absorption for the transitions of SiS (black), C2H2 (dark blue), H^^CCH (light blue) and HCN (green). 
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mation zone, placed at ~5 (jKeadv et al.lll993l ). We 
can derive the SiS abundance in that region from the 
v=l Hnes assuming an uniform sphere at a distance of 
180 pc with a radius of 10 R*, iUuminated by the central 
star (Rh,=5x1013 T,=2300 K), with Tk=1000 K and 
= 1.6x10^ cm~^ (mean valu es for this region derived 
bv lFonfria Exposito et al.]|2006f ) and an expansion veloc- 
ity of 11 km s""'^. We solve the statistical equilibrium for 
SiS considering 100 rotational levels and 3 vibrational 
states and apply the LVG radiative transfer formalism 
using the code developed by J. Cernicharo. With the as- 
sumed UHa, the column density for H2 is ~7xl0^'^ cm~^ 
and the derived one for SiS is ~5.0xl0^^ cm~^. Hence, 
the SiS abundance is ~7x 10^^. This result is compatible 
with a higher SiS abundance in t he innermost CSE (from 
LTE chemistry models, 3x 10~^. lAgundez fc Cernicharol 
I2OO6') and an abundance further away of 6.5x10""^ ac- 
cording to ob servations of v=0 rotationa l lines over the 
outer CSE bv lBieging fc Nguven-Q-Rieul (|1989D . 

Most v=0 rotational lines show a rounded or slightly 
double peaked profile with the blue part absorbed by 
cold SiS through the envelope. However, Mi, M2 and 
M3 show extra emission in the form of narrow peaks 
(FWHM=l-3 km s'^). The velocities of these fea- 
tures, within the —10 to 10 km s~^ range, indicate 
that SiS maser emission arises from different regions lo- 
cated between the star and the inner dust formation zone 
(r~5 R,.). The lines M 2 and M3 have been previously ob- 
served bv lSahai et al.l |i984) but no maser emission was 
noticed. This could be due either to the limited sensitiv- 
ity of their observations or to a time variability of the SiS 
maser phenomenon. The bottom-center panel of Fig. [1] 
shows in detail the line profiles of the observed masers. 
Up to ten maser features labelled (a,. . . ,j) are identified. 
The most complex line profile is that of M3. It is formed 
by 5 main features: a, d, /, i, j with v—S.S, —3.7, 
0.88, 6.3 and 8.4 km s~^, respectively. The line profiles 
show that the strongest peaks are at negative velocities, 
having their red counterpar ts rather weak . This behav- 
ior was previously found bv lHenkel et al.l (|1983D for the 
v=0 J=l-0 line, which mostly consists of a narrow peak 
(FWHM=0.3 km s"!) centered at -13.5 km s'^. The 
strong asymmetry of the lineshapes can be either due to 
blanking by the star of the redshifted maser or by am- 
plification of the blueshifted emission by the foreground 
stellar environment. 

The strongest observed maser, M3, with Tmb— 60 K 
(Fobs— 300 Jy; with thermal and non-thermal emis- 
sion), is weak comp ared to some SiO mas ers detected 
in 0-rich stars (e.g. ICernicharo et al.]|1993D or to HCN 
mase r s observed in IRC-l-10 216 ( Lucas fc Cernicharol 
[19891: iSchilke fc MentenI [20M ) . However, Mi, M2, and 
M3 are stronger than the SiS ■/=l-0 mase r observed to- 
wards IRC-H10216 bv lHenkel erall (|1983D . The similar- 
ity between maser features in M2 and M3 indicates that 
they may arise from the same regions and produced by 
the same pumping mechanism; the maser in Mi is prob- 
ably formed in other regions. Hence, we suggest two 
possible geometries of the innermost CSE to explain the 
observed features: 

(i) An onion-like innermost region, where each maser is 
produced in a shell. This hypothesis is supported by the 
symmetry of features a-j and d-i. The peaks at extreme 
velocities, a and j, would be produced just in front of 



and behind the star near the inner dust formation shell 
(r~5 R*) with expansion velocities of ~5-ll km s~^. The 
features d and i would be produced in a similar way but 
in an inner shell with a lower expansion velocity. Finally, 
the central peak, /, would be formed in a shell very close 
to the star, with the whole shell contributing to the maser 
emission. Mi would be produced in a cap-shaped region 
in front of the star. 

(m) All the masers are formed in different positions of a 
clumpy shell. The different features in M2 and M3 would 
be produced in different regions of the shell: peaks a and 
j in front of and behind the star and the other peaks 
(d, /, and i) in different clumps, as occur with the only 
feature of Mi . 

The classic pumping mechanism for the SiO v>0 
masers observed in 0-rich stars resides in the in- 
crease of the trapping lifetime (A/r)^_,„_i with J 
for f — »f-l transitions, when they become optically 
thick (|Kwan fc Scovillg |1974| ). Such mechanism pro- 
duces masers in adjacent rotational lines of the v 
state, and explains the v=l and 2 SiO 



masers 
Ehtzurj 



(i Buiarrabal fc Nguven-Q-Rieul Il981l : iLockett 
il992j ). However, the masers observed in rotational 
transitions of ^^SiO, ^°SiO, and in v=3 and 4 of SiO 
do not show the latter behavior and have been inter- 
preted as due to IR overlaps between ro-vibr a tiona l 
lines of SiO isotopologues (Cernicharo et al.l 119911 : 
iGonzalez- Alfonso fc Cernicharol 19971 ). For SiS. the ab- 
sence of maser emission in v=l rotational lines and the 
odd v=0 pattern also exclude the Kwan fc Scoville pump- 
ing mechanism. This suggests that overlaps of ?;=l-0 ro- 
vibrational transitions of SiS with those of mode 1/5 of 
C2H2 and mode 1^2 of HCN, could provide the pumping 
mechanism. C2H2 and HCN are abundant in the inner 
CSE of IRC+10216 and dominate the 11-14 /xm spec- 
trum (|Fonfrfa Exposito et al.ll2006h . Overlap s with these 
two sp ecies have been already proposed bv ISahai et al.l 
()1984f ) to explain the different profiles of adjacent J lines 
of SiS. However, the SiS frequencies used in that work 
were not as accurate (ct~10~^ cm~^) as those available 
today. We have calculated those fr equencies from th e 
Dunham coefficients determined by ISanz et all ([2003), 
for which the error of the band center is <10~^ cm~^ 
(~0.04 km s~^; the relative accuracy of P and R lines 
is much better). The frequencies of C2H2, H^'^CCH, 
HCN, and H^^CN fines have been taken fro m tfie HI- 
TRAN Database 2004 (jRothman et al.ll2005D . with an 
accuracy better than 10~^ cm"-'^ (~0.4 km s~^) for C2H2 
and H^CCH, and 10"'* cm"! for HCN and H^^CN. 

Table [T] shows the mid-IR line overlaps of SiS with 
C2H2, HCN, and their most abundant isotopologues. 
For the overlap search we selected coincidences within 
|Au|<10 km s^^. However, since the CSE is expanding, 
every region of the envelope is receding from the others. 
Hence, if the population of the SiS levels is affected by an 
overlap with a strong line of otfier species, the frequency 
of this overlapping transition must be higher than the 
SiS one. This would restrict the condition to positive 
Aw. Nevertheless, due to the linewidth, lines at At;<0 
can overlap the SiS lines. Hence, we fiave set the negative 
cutoff to one half o f the typical linewidth in the i nnermost 
CSE (=i5 km s"^: iFonfria Exposito et ani2006l ). There- 
fore, our search is restricted to —2.5 < A?; < 10 km s~^. 
All tfie ro-vibrational SiS lines commented fiereafter re- 
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TABLE 1 

Mid-infrared Line Overlaps of SiS with C2H2, HCN, and 
Their Most Abundant Isotopologues 



Line v (cm ^) 



Mol. 



Transition 



Ai'(km/s) 



Overlaps Involving SiS Levels of Observed Lines 



R 9 


750.3695 


HCN 


01^0 - 


- 00"0 Re (12) 


-6.6 


PIO 


738.2889 


C2H2 


1-Hi 


- 1-10° R/ (3) 


1.2 


R13 


752.6431 


C2H2 


0»li 


- 000° Re (9) 


5.9 


P14 


735.7324 


C2H2 


oOl-i 


- 0'^0° Qe (38) 


-6.1 


P15 


735.0861 


C2H2 




- 0<'00 Qe (36) 


-9.9 


Overlaps Invc 


Iving SiS Levels of Unobserved 


Lines 


P 2 


743.2624 


C2H2 


oOli 


- OOOO Re (5) 


0.6 


P16 


734.4368 


C2H2 


o"i-i 


- qOoO Qe (34) 


1.1 


P20 


731.8114 


C2H2 




- O^OO Qe (24) 


9.5 


P22 


730.4815 


C2H2 


l-Hi 


- 110° Qe (18) 


-2.1 


R22 


757.5819 


C2H2 


l^ll - 


- I^OO Re (10) 


0.1 


R22 


757.5819 


H"CN 


01^0 - 


- 00<'0 Re (17) 


5.2 


P23 


729.8123 


Hi^cCH O^l-i 


- O^OO Qe (19) 


8.7 


P24 


729.1402 


C2H2 




- OOOO Qe (1) 


9.4 


P25 


728.4654 


Hi^cCH 




- OOOO Qe (7) 


-0.0 


R25 


759.1733 


HCN 


01^0 - 


- 00"0 Re (15) 


3.6 


R26 


759.6977 


C2H2 


oOli - 


- 0"00 Re (12) 


-1.0 


R33 


763.2816 


h"cn 


01^0 - 


- 00<'0 Re (19) 


6.1 


P38 


719.4363 


C2H2 


l^li 


- l^OO Pe (5) 


4.5 


R40 


766.7130 


C2H2 


oOli - 


- O^OO Re (15) 


4.3 


R42 


767.6652 


C2H2 


1^1-1 


- l^O" Re (20) 


8.7 


R45 


769.0697 


C2H2 


oOli - 


- O^O" Re (16) 


-1.8 


P48 


712.1720 


C2H2 


1-H-i 


- l-iO° Pf (8) 


4.7 



Note. — Overlaps of SiS with C2H2, H"CCH, HCN and H^^CN 
found in the mid-IR range [690,780] cm^^ with -2.5 < An < 10 
km s"^ where Av / c = [i.'(X) - i/(SiS)] /i^(SiS) and J < 50, (corre- 
sponding SiS rotational frequencies below 900 GHz). The lines 
and frequencies at the left correspond to the vibrational transitions 
v—1 — >0 of SiS. The errors on the velocities are < 0.5 km s~^. The 
notation for the C2H2 and H^^CCH vibrational states involved in the 
ro-vibrational transitions is yt^^ , whereas for HCN and H^^CN is 
i/ii/|i/3. The parity of the lower level is even (e) and odd (/). The 
transitions are labelled as R, P, and Q for J^p — Jioxu — +1, —1, 0. 

fer to v—l^Q transitions and will be labelled with the 
usual spectroscopic nomenclature R, Q, P (see footnote 
of Tabled. 

In order to qualitatively interpret the effects of IR over- 
laps on maser emission, we have used the same LVG ra- 
diative transfer code modified to account for overlaps, 
changing the intensity at the overlapping frequency and 
the escape probability for photons from the overlapped 
lines. Thus, for some SiS lines, the excitation tempera- 
ture becomes negative (maser activation) and the bright- 
ness temperature is considerably enhanced. M2 and M3 



are naturally explained by the overlap of the R(13) line 
of SiS with the strong C2H2 fine OPl^ - 0°0°Re (9) (Ta- 
ble [T]). SiS molecules are easily excited from v—Q J=13 
to v=l J=14 through R(13) and decay to the v=Q J=15 
via P(15), creating a population inversion between the 
u=0 J=13 and 15, and producing maser emission in M2 
and M3. Ml may be produced by the overlap of the P(10) 
SiS line with the strong C2H2 line I'^l^ - I'^O^Ylf (3). 
This overlap can pump SiS molecules from v~Q J~IQ 
to v=\ J=9 through P(10), depopulates the f=0 J=10 
level and produces the inversion between v=Q J=10 and 
11. 

We have also looked for overlaps of u=l-0 higher- J SiS 
lines with C2H2 and HCN transitions to try to predict SiS 
masers at submillimeter wavelengths. Some of them are 
shown in the second block of Table [TJ They suggest, for 
example, that a maser could be found in rotational tran- 
sitions involving the v=Q J— 23 level (likely J=24-23 and 
maybe 23-22), or the i;=0 J=25 and 26 states (perhaps 
in the v=0 J=27-26 and maybe 26-25). These overlaps 
could also produce masers in v=l rotational transitions. 

There are many overlaps of other SiS isotopologues 
with lines of C2H2, HCN and ^^Si^^S. With the adopted 
criteria for the overlap search (see footnote of Table [1]) , 
we have found 91, 93, 76, and 94 coincidences for ^^SiS, 
^°SiS, Si^^'S, and Si^^S, respectively. Consequently, al- 
though these species are less abundant than ^^Si'^^S, the 
population of some levels could be inverted producing 
maser emission. 

This study represents the discovery of three new SiS 
masers and should be complemented with future observa- 
tions of higher- J v=0 and 1 rotational transitions. Fur- 
thermore, a detailed multi-molecule non-local radiative 
transfer model would help to understand the dust for- 
mation region and the role of SiS in C-rich evolved stars. 
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